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ABSTRACT 

Current models of galaxy evolution suggest that feedback from active galactic nuclei is needed to 
explain the high-luminosity cutoff in the galaxy luminosity function. Exactly how an AGN outflow 
couples with the ambient medium and suppresses star formation remains poorly understood. However, 
we have recently uncovered an important clue to how that coupling might work. Observations of Ha 
emission and blue light from the universe's most luminous galaxies, which occupy the centers of galaxy 
clusters, show that star formation happens only if the minimum specific entropy of the intracluster 
gas is < 30keVcm 2 . Here we suggest that this threshold for star formation is set by the physics of 
electron thermal conduction, implying that conduction is critical for channeling AGN energy input 
toward incipient star-forming regions and limiting the progress of star formation. 

Subject headings: 



1. INTRODUCTION 

What determines the upper limit to the luminosity 
of a galaxy? This question has received a great deal 
of attention lately because the answer seems to involve 
feedback from a galaxy's nucleus. Semi-analytic mod- 
els of galaxy formation show that feedback of some kind 
is necessary to account for the observed characteristics 
of the galaxy luminosity function. Models without feed- 
back overpredict the number of both low-luminosity and 
high- luminosity galaxies. Plausible amounts of feedback 
from supernovae bring predictions for low-luminosity 
galaxies into agreement with observations, but getting 
things right at the high- luminosity end is more difficult — 
the amount of supernova feedback needed to quench 
star formation in these systems is implausibly large 
(|Benson et al.l I2003D . Numerical simulations of galaxy 
formation treat the physics of merging and cooling more 
realistically but still do not solve the problem of the high- 
luminosity cutoff. In current simulations, late-time cool- 
ing of hot gas in the central galaxies of clusters produces 
so many stars at z ~ that the galaxies have luminosi- 
ties several times larger than the observed cutoff in the 
galaxy lu minosity function and colors that are much too 
blue fe.g.. lSaro etaHl2006t) . 

A closely related cooling problem has long bedeviled 
our understanding of galaxies and galaxy clusters, and 
X-ray observations strongly suggest that the solution lies 
with active galactic nuclei. In roughly half of nearby clus- 
ters, the cooling time at r < 100 kpc is less than a Hubble 
time, implying that gas at the center of the cluster must 
cool and condense, if there is no mechanism to offset cool- 
ing. However, the star- formation rates in these systems 
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are at most 1% to 10% of the expected gas condensation 
rate, and X-ray spectroscopy shows that the the bulk of 
the central gas remains suspen ded at temperatures j> 1/3 
of the virial temperature (see iPeterson fc Fabianl I2006L 
for a review) . High- resolution imaging with Chandra has 
revealed that many clusters and giant elliptical galax- 
ies with short central cooling times also have cavities in 
their ICM that coincide with radio-emitting plasma out- 
flows from the active nucleus in the central galaxy. In 
many cases the AGN energy output inferred from these 
cavities is similar to the X-ray cooling rate of the cen- 
tral gas — a strong piece of circumstantial evidence in fa- 
vor of AGN feedback as the mechanism th at limits cool- 
ing and star formation in high-m ass halos (jBirzan et al.l 
12004 iMcNamara k. Nulsenll2007f ). Not all clusters with 
short central cooling times have obvious cavities, but that 
finding can be understood if AGN heating is episodic. 

The most recent versions of semi-analytic models 
have therefore incorporated schematic implementations 
of AGN feedback and find that "radio-mode" feedback is 
a plausible explanation for both the observed cutoff in the 
galaxy luminosity function and the fact th at the largest 
galax i es in the universe are red, not blue (jCroton et al.1 
120061: iBower et al1l2006D . To match the observations, 
radio-mode feedback must occur preferentially in high- 
temperature systems and must provide enough feedback 
energy to suppress star formation without making many 
stars. In this context, AGN feedback is most effective in 
systems experiencing hot-mode accretion because an ex- 
tended hot galactic atmosphere is necessary to halt and 
thermalize the energy of the AGN outflow. 

Here we consider the implications of two recent studies 
showing that star formation in brightest cluster galaxies 
(BCGs) is closely linked to the entropy 1 structure of the 
intracluster medium (ICM). Section 2 presents evidence 
that star formation occurs only in BCGs whose central 
entropy is < 30keVcm 2 , equivalent to a cooling time 
< 5 x 10 s years. The star formation rates in these sys- 

As is customary in the X-ray cluster field, we quantify entropy 

using the adiabatic constant K = kTn e 2 ^ fsee lVofJ2005l . for more 
background) . 
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Fig. 1. — Dependence of Ha emission in BCGs on central entropy 
from Cavagnolo et al. (2008a). Central entropy Ko nas been mea- 
sured by fitting Chandra entropy profiles with a radial power law 
plus a constant, and Ljja values are a heterogeneous set taken from 
the literature. Orange filled circles show clusters with detected Ha, 
and black filled circles give upper limits for clusters without detec- 
tions (see Cavagnolo et al. 2008b for more details). Notice that 
strong Ha emission appears only when Ko < 30 keV cm2 ■ 

terns are generally much smaller than the gas cooling 
rates in the absence of feedback, and understanding why 
star formation still proceeds at a reduced rate is likely 
to tell us something about the coupling between AGN 
feedback and the ICM. In section 3 we suggest that the 
critical central entropy threshold for star formation is de- 
termined by electron thermal conduction, because multi- 
phase structure cannot persist in the ICM if the central 
entropy is much larger than the observed threshold. If 
conduction is indeed responsible for governing the rate 
of star formation, then it is also likely to be a major con- 
duit of thermal energy in the cluster core and a crucial 
part of the AGN feedback loop. Section 4 summarizes 
the paper. 

2. THE STAR-FORMATION THRESHOLD IN BRIGHTEST 
CLUSTER GALAXIES 

Star formation in brightest cluster galaxies has long 
been known to correlat e with the centr al cooling time 
of the host cluster. iHu et all |l985) showed that 
Ha nebulosity is present in a BCG only if the clus- 
ter's cooling time is less than a Hubble time, and 
much of the Ha emission seen in these objects can 
be attributed to pho t oionization by hot stars (e.g. , 
Johnstone et all 119871 : iMcNamara fc O'ConneU 119891 : 
Voit &: Donahue! fl997f h While it is possible that some 
of the Ha comes from other sourc es like conduc- 
tive interfa ces (e.g., [Spa rks et al. 2004]), turbulent mix- 
ing layers (iBegelman fc Fabianl 119901) or co smic rays 
(|Ruszkowski et all 120081: iFerland et al l 12008ft . there is 
a clear connection between Ha emission and the pres- 
ence of extended ex cess blue light within the BCG (e.g., 
iCardiel et al.| [l998). Furthermore, clusters with short 
central cooling times also often have substantial CO 
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Fig . 2. — Depend e nce o f BCG color gradient on central entropy 
from [Raffertv et al. (2008). Color gradients have been measured 
in U — I and U — R, with open circles representing G(U — I) = 
d(U — I)/dlogr and filled circles representing G(U — R) = d(U — 
R) / d\ogr . The dashed line indicates a zero color gradient. The 
only BCGs with unusually blue centers are those in clusters with 
central entropy < 20 keV cm 2 , indicating that star formation is 
present only in those BCGs below this central entropy threshold. 

luminosities, indicating that they contain enough cool 
molecular gas to sup port the observed level of star- 
formation (jEdgdl2001l ). 

Two recent studies have brought the star-formation 
threshold into sharper focu s. One is a Chandra arc hival 
survey of galaxy clusters by I Cavagnolo et all (|2008l ) that 
measured the radial entropy profiles of 222 galaxy clus- 
ters, fitting them with the three-parameter expression 



K(r) = K + K 



100 
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Here we will refer to Kq as the central entropy, although 
in some cases it is not necessarily the minimum entropy 
level of the ICM but rather the entropy scale at which 
the radial profile departs from the best-fitting power law 
at large radii. Figure [1] shows that central Ha emission 
has been detected only in clusters with K q < 30 keV cm 2 
(Cavagnolo et al. 2008b). Emission- line luminosities in 
this plot have been taken from the literature and corre- 
spond to a variety of apertures, so the quantitative values 
of Ln a are not very meaningful. However, the entropy 
threshold for having detectable Ha emission is quite dis- 
tinct. The presence of Ha emission clearly indicates that 
the intracluster medium can have a multiphase struc- 
ture below this entropy threshold, with a component at 
~ 10 4 K in addition to the X-ray emitting component at 

> 10 7 K. 

Observations of color gradients in BCGs have led 
to a similar conclu sion. In a sample of 46 clusters, 
iRaffertv etall (|2008l ) find that 34 of the 20 clusters with 
central entropy < 30keVcm 2 have BCGs with extended 
blue light in their cores consistent with ongoing star for- 
mation, while none of the clusters with central entropy 

> 30keVcm 2 have BCGs with colors bluer than expected 
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from a passively evolving stellar population (see Fig- 
ure [2]) Here the central entropy has been measured from 
a small X-ray aperture ~ 2 — 3 arcsec at the center of the 
cluster, and the color gradients have been measured in 
U - 1 and U - R, with G{U - I) = d(U - I) jd log r and 
G(U -R) = d(U - R)/d\ogr. 

The transition to Ha emission and star formation is 
sharp in both cases, but it is not yet clear whether the 
threshold ought to be expressed in terms of central en- 
tropy or central cooling time. Here we are presenting 
the results as a function of central entropy because we 
are going to interpret them in terms of the ICM entropy 
structure in the next section. However, it is interesting 
that the central cooling times at which Ha and blue light 
appear are substantially smaller than 10 10 years. Sim- 
ply having gas that is able to cool within a Hubble time 
is not enough for star formation in a BCG. Some other 
more restrictive condition must also be satisfied. 

It is also interesting that in both figures there are clus- 
ters with central entropy < 20 keV cm 2 and cooling time 
< 3 x 10 8 years that do not show evidence for star forma- 
tion. Low entropy and a correspondingly short cooling 
time are therefore not sufficient to guarantee star forma- 
tion. Rafferty et al. (2008) present evidence suggesting 
that star formation proceeds only if the AGN heating 
rate is less than the gas cooling rate and the central 
galaxy is close to the X-ray brightness peak. 

3. CONDUCTION AND THE CONDITIONS FOR A 
MULTIPHASE ICM 

We suspect that the critical entropy threshold for mul- 
tiphase gas and star formation in BCGs may result from 
electron thermal conduction. Cool star-forming clouds 
should appear only in systems whose size is greater than 
a critical length scale, known as the Field length, be- 
low which therm al conduction smoothes out temperatur e 
inhomogeneities (|Fieldlll965t iBegelman fc McKedll990l) . 
One can derive the Field length heuristically by consider- 
ing thermal balance for a cool cloud of radius r embedded 
in a medium of temperature T. Electron thermal conduc- 
tion sends energy into the cloud at a rate ~ r 2 n(T) ■ T/r, 
where k(T) = 6 x 10" 7 T 5 / 2 / c ergs s^K^cm -1 is the 
Spitzer conduction coefficient and f c is a suppression 
factor depending on the magnetic field structure in the 
medium. Radiative cooling can rid the cloud of en- 
ergy at a rate ~ r 3 n 2 A(T), where the cooling function 
A(T) cx T 1 ' 2 for T > 2 keV. Cooling and conduction 
are therefore in approximate balance for systems with a 
radius of order the Field length, 
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Through a coincidence of scaling, the Field length is a 
function of entropy alone whe n free-free emission i s the 
dominant cooling mechanism (|Donahue et al.ll2005l ). 

Figure 3 illustrates how this criterion translates into 
the entropy-radius plane. The long-dashed lines give the 
loci of points for which Xp(K) = r, given suppression 
factors f c — 0.2 and 1. Magnetic suppression of conduc- 
tion is a complicated and incompletely understood pro- 
cess, but most recent estimates h ave been in the range 
f c » 0.2-0.3 (jMalvshkinl 120011: iNaravan fc Medvedevi 
l2001h . Below each line, gas within radius r constitutes 
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Fig. 3. — Regions of thermal stability in the ICM based 
on the Field-length criterion. Long-dashed lines show where 
Xp(K) = r for conduction suppression factors f c = 0.2 and 1. 
Above these lines, conduction should be more effective than ra- 
diative cooling, causing cooler structures of scale < r to evap- 
orate. Below these lines, radiative cooling should be more ef- 
fective than conduction, allowing thermal instability to proceed. 
Solid lines show schematic entropy profiles of the form K(r) = 
K + (150keVcm 2 )(r/100kpc) 1;2 , with K = 0, 10, 20, 30, and 
50, as labeled. Short-dashed blue lines give entropy profiles from 
(Cavagnolo et al. 2008a) for clusters in which Rafferty et al. (2008) 
find central star formation. Dotted red lines give entropy profiles 
for clusters from Rafferty et al. (2008) without clear evidence for 
star formation and that also have no detectable Ha emission (Cav- 
agnolo et al. 2008b). 

a subsystem with r > Ap. Conduction cannot stabilize 
that gas against cooling, allowing multiphase gas to per- 
sist and star formation to proceed. Above each line is the 
region of stability, in which conduction leads to evapora- 
tion and homogeneity. 

For comparison, the solid lines in Figure 3 show 
schematic intracluster entropy profiles motivated by 
the large Chandra entropy survey of Cavagnolo et al. 
(2008a). Fitting the form in equation to those profiles 
gives K 100 = 146.8 ± 68.1 keVcm 2 and a = 1.22 ± 0.26 
for clusters with Kq < 50keVcm 2 and temperatures 
ranging from 2 keV to 12 keV. The solid lines there- 
fore represent typical intracluster entropy profiles with 
K wo = 150keVcm 2 , a = 1.2, and K = 0, 10, 20, 30, 
and 50 keV cm 2 . Notice that only the schematic profiles 
with Kq < 30keVcm 2 , corresponding to clusters with 
central star formation and Ha emission, dip below the 
threshold for conductive stabilization corresponding to 
f c — 0.2. Colored lines giving entropy profiles from (Cav- 
agnolo et al. 2008a) reinforce this correspondence: short- 
dashed blue lines show Rafferty et al. (2008) clusters 
with central star formation, and dotted red lines show 
Rafferty et al. (2008) clusters without clear evidence 
for star formation and without detectable Ha emission 
(Cavagnolo et al. 2008b). 

This simple analysis suggests that conduction is re- 
sponsible for the strong dependence of Ha emission and 
star formation on central entropy. Even though the sup- 
pression factor f c is uncertain, a threshold must exist in 
the K-r plane above which conduction heats and evap- 
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Fig. 4. — Dependence of BCG color gradient on conduction 
efficiency. Open and filled circles have the same meanings as in 
Figure [5] Conduction efficiency is quantified by finding the mini- 
mum value of kT / (A(T)n e nnr 2 ) for each cluster. If this quantity 
is less than l// c , then conduction cannot offset radiative cool- 
ing. The fact that blue gradients appear only in clusters with 
kTJ (A(T)n e nnr 2 ) < 5 somewhere in the cluster suggests that ther- 
mal conduction with f c ss 0.2 determines whether a long-lasting 
multiphase medium can be present in the cluster core. 

orates cooler gas clouds faster than radiative losses can 
cool them. One therefore expects to see a transition from 
a multiphase ICM at low central entropy levels to a more 
homogeneous, isothermal ICM at higher entropy levels. 
The Ha observations are indicating that this transition is 
at approximately the expected location in the K-r plane 
for thermal conduction with f c w 0.2. We note that the 
presence of such a threshold is independent of the origin 
of the line-emitting gas, which could come either from 
condensation of the ICM or from stripping of cool inter- 
stellar clouds from galaxies passing through the cluster 
core. Stripped gas in an ambient medium above the crit- 
ical K(r) profile would quickly evaporate, while stripped 
gas in a medium below the profile could persist indefi- 
nitely. 

A case-by-case analysis of the effectiveness of conduc- 
tion in the clusters from iRaffertv et al.l (|2008f ) agrees 



with the broader conclusions drawn from Figure 3. For 
each cluster we have determined the minimum value of 
KT/[A(T)n e riB_r 2 }, where A(T) represents the full cooling 
function and not just free-free emission, from Chandra 
observations. If KT/[A(T)n e n-ar 2 ] < f^ 1 , then conduc- 
tion cannot compensate for radiative cooling (see Figure 
4). The fact that blue gradients are seen only in the 
BCGs of clusters with KT/[A(T)7i e n H r 2 ] < 5 therefore 
provides further support for the idea that thermal con- 
duction with f c w 0.2 determines whether or not a BCG 
can form stars. 



4. SUMMARY 

Two recent studies have shown that the star-forming 
properties of BCGs depend on the entropy structure 
of the host cluster's ICM. Clusters with central en- 
tropy < 30 keV cm 2 , corresponding to a cooling time 
< 5 x 10 8 years often have BCGs with Ha emission 
(Cavagnolo et al. 2008b) and blue cores (Rafferty et al. 
2008) indicative of ongoing star formation, while none of 
the clusters with high central entropy show the same fea- 
tures. We suggest that this entropy threshold for star for- 
mation is set by the physics of thermal conduction. Clus- 
ters typically have entropy profiles that approximately 
follow the law K (r) = K + (150keVcm 2 )(r/100kpc) 1 ' 2 , 
with differing values of Kq (Cavagnolo et al. 2008a). If 
Kq > 30keVcm 2 , then conduction can carry heat en- 
ergy toward the center of the cluster faster than the 
ICM can radiate it away, as long as the suppression fac- 
tor fc <; 0.2. Therefore, we should not expect clusters 
with higher central entropy levels to have a persistent 
multiphase medium in their cores. Conversely, conduc- 
tion in lower-entropy clusters cannot fully compensate 
for radiative coo ling, allowing star formation to proceed. 
IGuo et aH (|2008t ) have recently presented a model ex- 
plaining how AGN feedback can suppress star formation 
in these thermally unstable cluster cores. The corre- 
spondence we find between the observed star-formation 
threshold and the predictions of conductive models in- 
dicates that conduction indeed is important in cluster 
cores and may be crucial to the AGN feedback mecha- 
nism thought to regulate star formation in these systems. 
However, detailed hydrodynamical simulations including 
conduction will be necessary to determine more precisely 
how the phase structure of the ICM depends on a clus- 
ter's entropy profile. 
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